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Abstract 

The statistical distribution of eigenvalues of pairs of coupled random matrices can be expressed in terms 
of integral kernels having a generalized Christoffel-Darboux form constructed from sequences of biorthogonal 
polynomials. For measures involving exponentials of a pair of polynomials Vi , V2 in two different variables, these 
kernels may be expressed in terms of finite dimensional "windows" spanned by finite subsequences having length 
equal to the degree of one or the other of the polynomials Vi, V2. The vectors formed by such subsequences 
satisfy "dual pairs" of first order systems of linear differential equations with polynomial coefficients, having 
rank equal to one of the degrees of Vi or V2 and degree equal to the other. They also satisfy recursion relations 
connecting the consecutive windows, and deformation equations, determining how they change under variations in 
the coefficients of the polynomials Vi and V2. Viewed as overdetermined systems of linear difference-differential- 
deformation equations, these are shown to be compatible, and hence to admit simultaneous fundamental systems 
of solutions. The main result is the demonstration of a spectral duality property; namely, that the spectral curves 
defined by the characteristic equations of the pair of matrices defining the dual differential systems are equal upon 
interchange of eigenvalue and polynomial parameters. 
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1 Introduction 



1.1 Random matrices 

1.1.1 Background and motivation 

Random matrices |]37| , Q play an important role in many areas of physics. They were first introduced 
by Wigner [|T) in the context of the spectra of large nuclei, and the theory was greatly developed in 
pioneering work of Mehta, Gaudin and Dyson |l^. It has found subsequent applications 

in solid state physics [p^ (e.g., conduction in mesoscopic devices, quantum chaos and, lately, crystal 
growthp^), in particle physics |5^, 2d-quantum gravity and string theory jll], [ij, ||. The reason for the 
success and large range of applications of random matrices is, due, in part, to their universality property; 
when the size of the matrices N becomes large, the statistics of the eigenvalues tend to be independent of 
the model, and determined only by its symmetries and the spectral region considered, relative to critical 
points and edges in the spectral density. Matrix integrals are also known to give special realizations 
of KP, Toda and isomonodromic r-functions, and thus have a close relationship to integrable systems 
@,|6|,||,0,|,|,||,§. 

Random matrices also have important applications in pure mathematics, for example, in the statistical 
distribution of the zeros of functions They are also related to other statistical problems 

such as random word growth and the lengths of nondecreasing subequences of random sequences Q . 

The model we shall consider here is called the "2-matrix model" |8[ [l^, This 
involves an ensemble consisting of pairs oi N x N hermitian matrices Mi and M2, with a U{N) invariant 
probability measure of the form: 

—dn{Mi,M2) := — exptr (-V^i(Mi) - V2{M2) + MxM2)dMidM2 , (1-1) 

where dMidM2 is the standard Lebesgue measure for pairs of Hermitian matrices, Vi and V2 are polyno- 
mials of degrees di -I- 1, d2 + l respectively, called the potentials, with coeffficients viewed as deformation 
parameters, and the normalization factor (partition function) is 



TN= dlJi, (1-2) 

J Ml JM2 

which is known to be a KP r-function in each set of deformation parameters, as well as providing solutions 
to the two- Toda equations ^, |^. 

This model was introduced in ||30|, p8| as a toy model for quantum gravity and string theory. The 
main interest was in a special "double scaling" limit, where N ^ 00 and the potentials Vi and V2 are 
fine-tuned to critical potentials. The asymptotic behaviour in such limits is related to finite dimensional 
irreducible representations of the 2D-conformal group |lj, |lo) . The best known example is when Vi 
and V2 are cubic polynomials, tuned to their critical values, which reproduces the critical behaviour of the 
Ising model on a random surface ||3^, js) . It is important to note that the 2-niatrix model contains more 
critical points than a 1-matrix model; for instance, the 1-matrix-model cannot have an Ising transition 

0- . . . ^ ^ 

String theorists have also introduced a generalization, known as the "multi-matrix model" ||l3|, y^, 

H, m, where one has a set of m > 2 matrices {N x N hermitian) coupled together in a chain, with a 
measure of the form 

—dfi{Mi , . . . M„) = — exp tr { ~^Vj (Mj) + ^ M,M,+i [| dM, , (1-3) 



tn tn 
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and the V^ 's are again polynomials in their arguments. This model has the same universal behaviour as 
the 2-matrix model and, in some sense, does not seem to contain any more information. Throughout 
the main body of this work, we will concentrate on the 2-matrix-model, for which the statistics of the 
eigenvalues can be calculated using biorthogonal polynomials ^ |2^, ||, |^ . In the appendix, it will 
be explained how to extend all the results in the present work from the 2-matrix model to the differential 
systems associated with this multi-matrix model. 



1.1.2 Relation to biorthogonal polynomials 

By biorthogonal polynomials, we mean two sequences of monic polynomials 

7r„(x) =a;"-f ct„(j/) = y" + • • • , n = 0, 1, . . . (1-4) 

which are orthogonal with respect to a coupled measure on the product space: 

dxdy ^„(:E)a,„(y)c-^i(-)-^^(^')+^y = /i„5,^„, (1-5) 



where Vi ( x) a nd V2 (y) are polynomials chosen to be the same as those appearing in the 2-matrix model 
measure (1-1), and a suitable contour is chosen to make the integrals convergent. The orthogonality 
relations determine the two families. Once the biorthogonal polynomials are known, they may be used 
to compute four different kernels: 



Ki2ix,y) = J2 7^7r„(a;)a„(y)e-^i(-)e-^^(^) , Kii{x,x') = / dy Ki2ix,y) e^'y , (1-6) 

n— 

N f N , N f f N ' ' 

K22{y\y) = dx Ki2{x,y)e''y , K2i{y,x')^ dxdyKi2{x,y) e^^^ e"" ^ . (1-7) 

All the statistical properties of the spectra of the 2-matrix ensemble may then be expressed in terms 

mi ^ 
of these kernels |22] and the corresponding Fredholm integral operators Kij, i,j = 1,2. For instance the 

density of eigenvalues of the first matrix is: 

N IN 

Plix) ^ —Kii{x,x) , (1-8) 
the correlation function of two eigenvalues of the first matrix is: 

N , I / N N , , N , ^ , \ 

Pn{x,x') = ^lKii{x,x)Kii{x',x')-Kii{x,x')Kii{x',x)] , (1-9) 
and the correlation function of two eigenvalues, one of the first matrix and one of the second is: 

N I / N N N N \ 

Pi2{x,y) ^ — iKiiix,x)K22{y,y) - Ki2ix,y)iK2i{y,x) - e^y)\ . (1-10) 

Any other correlation function of m eigenvalues can similarly be written as a determinant involving these 
four kernels only. 
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The spacing distributions (the probabihty that two neighbouring eigenvalues are at some given dis- 
tance) can be computed as Fredholm determinants. For example, the probability that some subset J of 
the real axis contains no eigenvalue of the first matrix is the Fredholm determinant: 

p^-' = det (^1 - Kii o xjj , (1-11) 

where XJ is the characteristic function of the set J. 

An important feature in the study of the N oo limit is that the kernels Kij may be expressed [^o| 
in terms of sums involving only a fixed number of terms (either di + 1 or ^2 + 1), independently of N, 
as a consequence of a "generalized Christoffel-Darboux" formula ^ following from the recursion 
relations satisfied by the biorthogonal polynomials. This allows one, in the N ^ oo limit, with suitable 
scaling in the spectral variables, depending on the region considered, to treat TV as just a parameter. 



1.2 Duality 

1.2.1 Dual isomonodromic deformations 

The notion of duality arises in a number of contexts, both in relation to isospectral flows [Q and isomon- 
odromic ||2^, ^ deformations. What is meant here by "duality" in the case of isomonodromic defor- 
mations is the existence of a pair of parametric families of meromorphic covariant derivative operators 
on the Riemann sphere 

V^:^ — +L{x,u), I?2 := ^ + M(y,u) , (1-12) 

where L(a;,u) and M(j/,u) are, respectively, / x I and m x m matrices that are rational functions of the 
complex variables a; € and y € P^, with pole divisors of fixed degrees, depending smoothly on a set 
of deformation parameters u — (ui, U2, ■ ■ ■) in such a way that: 

1) The matrices L(a;,u) and M(?/, u) are obtained from the integral curves of a set of commuting (in 
general, nonautonomous) vector fields defined on a phase space M by composition with a prescribed pair 
of maps (possibly depending explicitly on the deformation parameters) from A4 to the spaces of rational, 
I X I or m X m matrix valued rational functions of the spectral parameter x or y, respectively, with pole 
divisors of fixed degree. 

2) The generalized monodromy data of both the operators I?i and I?2 are invariant under the u- 
deformations. (This includes the monodromy representation of the fundamental group of the punctured 
Riemann sphere obtained by removing the locus of poles and, in the case of non-Fuchsian systems, the 
Stokes matrices and connection matrices ^.) 

3) The spectral curves determined by the characteristic equations: 

det(L(a;,u) -2/1) = , det(M(y, u) - xl) = (1-13) 
are biholomorphically equivalent. 

(A similar definition can be given for the case of dual isospectral fiows of matrices L(a;, u) and M(2/, u).) 

Such "dual pairs" of isomonodromic families occur in many applications. They are related to the 
solution of "dual" matrix Riemann-Hilbcrt (RH) problems ^ ^ which, in certain cases, are equiv- 
alent to determining the resolvents of a special class of "integrable" Fredholm integral operators K, K, 
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with kernels of the form: 



^^MEkiM, (1-14) 

^ — ^ X — X 

i—1 

t{ y-y' 



Here, the vector valued functions 



= (/i(x,u),.../;(a;,u)) ; g(a;,u) = {gi{x,u), . . . gi{x,u)) (1-16) 
f(?/,u) = (/i(?;,u),.../„(y,u)) ; g(y, u) = (gi(y, u), . . . 5,„(?/, u)) . (1-17) 

depend on the spectral variables x and y as well as on some, but not necessarily all the deformation 
parameters {ui,U2, ■ ■ ■)■ They satisfy overdetermined, compatible differential systems in these variables 
which imply the invariance of the monodromy of associated "vacuum" isomonodromic families of covariant 
derivative operators 

25oa := ^+Lo(a;,u) , I?o,2 ^ + Mo(y, u) . (1-18) 

The dual kernels K{x,x') and K{y,y') are related to each other by applying partial integral transforms 
with respect to one of the two spectral variables {x,y) (e.g., Fourier-Laplace transforms) to an integral 
operator on the product space x P^. Application of the Riemann-Hilbert dressing method for suitably 
chosen "dual" sets of discontinuity data then gives rise to the "dressed" families ( 1-12[ ), which have 



a similar relation to the resolvent kernels of the two operators. (Additional deformation parameter 
dependence may enter, besides that contained in the vacuum equations, characterizing the support of 

these operators.) The Fredholm determinants det(l — K) and det(l — K) may then be shown, through 
deformation formulae, to coincide with the corresponding isomonodromic tau functions, and with each 
other (see e.g. ||, H]). 

Such systems arise naturally, as discussed above, in the study of the spectral statistics of random 
matrix ensembles, both in the finite case, and in suitable infinite limits. An example of such dual pairs 
is given by the class of exponential kernels in which: 

Mx) - (-l)'e"-- , 5,(x) = e-"'- , (1-19) 
fa{x) - i-ire^^- , U^) = e-""- , (1-20) 

For the case I = m — 2, these include the sine kernel 

Ksix,x')^'J^p^, Ksiy,y')^'^^f^ (1-21) 
(X - x') (y - y') 

governing the spectral statistics in the scaling limit of the GUE in the bulk regio n |37|. 

Other examples include the various Painleve equations P//, P/y, Py, Pvi |26|, 12^ which each 
possess "dual" isomonodromic representations. A last class of examples, unrelated to random matrices, 
but including a special case of Py/, consists of the isomonodromic representations of the WDVV equations 
of topological 2D gravity entering in the theory of Frobenius manifolds |l5) . These possess both Fuchsian 
representations with n finite poles with residues of rank 1, and non- Fuchsian n x n representations, having 
a single irregular singular point of Poincare index 1 at cx). 
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1.2.2 Duality in the large N Limit 

It is proved in under a suitable large N assumption, and conjectured for other cases ^ 

that in a particular large N limit (where the coefficients {uk, vj} scale as N, and the supports of 

N N 

Pi{x) :— lim7v->ooPi(2^) and p2{y) '■= lim^v -^00/^2(2/) are connected intervals [ai,6i], [02,52] respectively), 
the following functions: 

/(.) := lv,'(.) - /" . - £• 01ay' , (1-22, 



are inverses of each other. In other words, if y = f{x) then g{y) — x. One can see that the functions f{x) 
and g{y) are related to the eigenvalues of the operators which implement the derivative with respect to 



X and y for the biorthogonal polynomials. The spectral duality theorems 4.1 and 4.2 which we present in 



this work are a more precise statement related to this conjecture, with a rigorous proof valid for all N . 
1.3 Outline of the article 

1.3.1 Biorthogonal polynomials and differential systems 

In Section 2, we consider the normalized quasi-polynomials 

V'„(:r) = -^7r„(a;)e-^^(^) , <^„(y) = ^a„(2/)e-^^(^) , n = 0, . . . oo . (1-23) 



Viewing these as the components of a pair of column vectors 

4- = (V'o,'(/'i,---V'n,---)* and $ = (00,01, •)* > (1-24) 

OO OO 

we obtain a pair of semi- infinite matrices Q and P that implement multiplication of ^ by x and derivation 

OO 

with respect to x, respectively. Equivalently, we obtain the transposes and by applying — ^ or 
multiplication by — y to By construction, these satisfy the Heisenberg commutation relations 

OO 

[P, Q] = 1 , (1-25) 

and, as shown in Section 2, they are finite band matrices; Q has nonvanishing elements only along 
diagonals that range from 1 above the principal diagonal to d2 below it, and P has nonvanishing elements 
only along the diagonals from 1 below the principal to c?i above it, where di -I- 1 and d2 + l are the degrees 
of the polynomials Vi{x) and V2(y), respectively. 



The first result (Prop. 2.1) following from the finite recursion relations satisfied by the quasi-polynomials 
{'>pn{x)} and {(f)n{y)} is a set of "generalized Christoffel-Darboux relations |2^, which imply that the 

N , ^ , 

kernels Kii{x,x ) and K22{y ,y) may be expressed as: 

/N-l N \ fN-l N 

^ t {x'),A^{x)] ^ * (?/'), B$(y) 
Kiiix, x') = - A , K22{y', y) = ^ ; , (1-26) 

X — x' y y 

where vl/(x) and $(2/) are the ^2+1 and di+1 dimensional column vectors with components [ipN~d2 , ■ • ■ , V'tv] 

N N 

N-l N-l 

and [<j)N~d2 , • ■ • , 0Af] , respectively, and ^ (y) and $ (x) are the d2 + l and di + 1 component row vectors 
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with components • ■ • , and [?!'^_-^, ■ • ■ , '^j^^^ respectively, where the underbarred quan- 

tities {'0 (y)} and {(f> (y)} designate the Fourier-Laplace transforms of the quasi-polynomials {ipn{y)} 

N N 

and {(/)„(?/)}. The matrices A and B are, essentially, the nonvanishing parts of the matrices obtained by 
commuting Q and P, respectively, with the projectors to the appropriate finite-dimensional subspace. A 
similar "differential" form of the generalized Christoffel-Darboux relations holds, following from applying 

N ^ , 

the derivations dx + dx' and dy + dy' to the kernels Kii{x, x ) and K22{y , v)- 

The recursion relations satisf ied by the quasi-polynomials {ipn{x)\n=Q...oo and {4>n{y)}n=o.. .00 may be 
conveniently expressed (Lemma 2.3) as: 



NN 



* {x) 

N+l 



b$(2/) 

N N 



$ {y) , 

JV+l 



where the "ladder" matrices a and b are linear in x and in w, respectively, and are formed from the rows 

N N 



of Q and the columns of P (see eqs. ( 2-4g| , 2-50 ) for their exact definitions) 



The vectors "^{x) and $(j/) also satisfy the following differential equations (Lemma p. 4| 

N N 



d 

— * 

dx N 



N 

Dl{x)^ll 

N 



d 

— $ 

dy N 



N 

D2{y)^ 

N 



(1-27) 



N N 

where the matrices Di{x) and D2{y) are, respectively, of size {d2 -f 1) x (^2 + 1) and {di -I- 1) x {di + 1), 
with entries that are polynomials in the indicated variables of degrees di and c?2 , respectively (and also 
polynomials in the matrix entries of Q and P). Furthermore, if {uK}K=i...di+i and {wj} j=i...d2+i are 
the coefficients of the polynomials V{{x) and V^iy), respectively, and these are varied smoothly, the effect 
of such deformations is given by the following system of PDE's (Lemma 2.S) 



d 



duK 
d 



N 



duK N 



TV,* 

U X* 

N 

U x$ 

N 



d 

V]/ 

dvj N 
d 

$ : 

dvj N 



- V j^- 

N 

N.'S' 

V J $ , 

N 



(1-28) 



where the matrices U xix), V j{x), U xiy) and V j{y) are again polynomials in the indicated vari- 
ables and in the matrix entries of Q and P. 



1.3.2 Compatibility 

So far, these statements are just a re-writing of the infinite series of recursion relations, differential 
equations and deformation equations satisfied by the functions {ipn{x)}n=o...oo and {<t>niy)}n=o...oo, pro- 
jected onto the finite "windows" represented by the vectors \1/ and However, we may now view these 

N N 

equations as defining an overdetermined system of finite difference-differential-deformation equations for 
vector functions of the variable {TV, x, y, Ui^, f j}, and ask whether, as such, these systems are compat- 
ible; i.e., whether they admit a basis of simultaneous linearly independent solutions. The affirmative 



answer to this question is provided in Section 3 by Prop. 3.3, which states that sequences of invertible 
{d2 + 1) X {d2 + 1) and (di -I- 1) x (di + l) matrices ^(x) and $(y) exist (fundamental solutions), for which 

N N 

all the column vectors satisfy the above difference-differential-deformation equations simultaneously. 

The compatibility of the deformation equations and finite difference equations with the x and y 
differential equations imply, in particular, that the (generalized) monodromy of the polynomial covariant 
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derivative operators 



ON d ^ 



(1-29) 



is invariant under both the {uktVj} deformations and the shifts in N. A similar statement can be made 
of the corresponding operators 

d ^ , . d 



N 



dx 



dy 



DM 



(1-30) 



N-l N-1 

annihilating the Fourier-Laplace transformed vectors $ (a;) and (y). 
1.3.3 Spectral duality 



A result that is not at all obvious from the above discussion is proved in Prop. 4.1; namely, that the 

N N N N 

pairs of matrices {Di{x), D.2{y)) and {D2{y) , D.iix)) have the same spectral curves. More specifically, we 
have the following equalities between their characteristic equations, which actually are identities 



Mdi+i det 
Mrfi+i det 



N 

xid,+i - D2{y) 



N 

DM 



= Vd^i det 
= Vd.^i det 



N 



yld^+i-Di{x) 



N 



yld2+i - Di{x) 



(1-31) 
(1-32) 



(Note that the two integers di + I and d2 + 1 which determine, in one case, the dimension of the matrix, 
in the other, its degree as a polynomial in the variable x or y, are interchanged in these equalities, as are 
the roles of the variables x and y.) 

What is even less obvious, and actually depends on the validity of the Heisenberg commutation relation 
satisfied by Q and P, is that these curves are not only pairwise equal, but in fact they all coincide, since 

N N N N 

the pairs of matrices {Di{x)), D.i{x)) and (£'2(2/)); Djziv)) ^'"^ conjugate to each other (Theorem | 
so we also have the equalities: 



det 
det 



N 

-Di(x) 

N 

xid^+i - D2{y) 



= det 
= det 



2+1 



N 

D,{x) 

N 

xld,+i -D^iy) 



(1-33) 
(1-34) 



Moreover, the transformations relating them are x and y independent, and are just those defined by the 
matrices entering in the generalized Christoffel-Darboux relations: 



N N N N 



N N N N 

MD2iy) = D2{y)M 



(1-35) 



These same matrices also relate the system of deformation equations, where they enter as gauge trans- 
formations depending on the deformation parameters {uk,vj}. 

The key to proving all these results lies in noting (Theorem |4.lD that, as a consequence of the differen- 
tial equations and recursion relations satisfied by the ipn^s, (pn^ and their Fourier-Laplace transforms, the 
following quantities are in fact independent of all the variables N , x, y, {uK}K=i...di+i and {vj}j=i,,,d2+i- 



'N-l 



N ~ 

* (j/),B|(y) 



9n{x) 



N ~ 

$ (a;), A*(a;) 



(1-36) 
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iV-1 Af-1 

where ^{x), ^(y), ^ {y), ^ (a;) are any solutions of the full system of difference-differential-deformation 
equations. This allows us to conclude that there exist compatible sequences of fundamental solutions \I>(x) 

N 

and $(y) of the above difference-differential-deformation equations, and the corresponding equations for 

N 

the Fourier-Laplace transformed quantities ^ , ^ such that 

JV-l AT \ /N-l N \ 

*,A*)j=l, f *,B*j=l, (1-37) 

for all values of {N, x, y, uj, vk}- 

This fact may be viewed as a form of the bilinear identities implying the existence of t- functions ^ . 
The development of this relation, and its connection with the isomonodromic deformation equations, 
which requires a study of the formal asymptotics of the fundamental solutions near a; = oo and y = oo 
will be left to a later work . 

In the appendix, all the above results are extended to the sequences of multiorthogonal functions that 
replace the biorthogonal quasi-polynomials in the multi-dimensional case associated to the multi-matrix 
model discussed above. 



2 Biorthogonal Polynomials 
2.1 Biorthogonality measure 

We first consider sequences of biorthogonal polynomials with respect to the measure arising in the study 
of the two-matrix model discussed in the introduction. Consideration of the recursion relations obtained 
by multiplication by or derivation with respect to the independent variable gives rise to representations 
of the Heisenberg relations (string equation) in terms of pairs of semi-infinite matrices. However, most 
of the results obtained here may also be shown valid in the fully infinite case. 
Let us fix two polynomials, which we refer to as the "potentials" , 

^i(-)=E^-^' ^2(y)=EV- (2-1) 

K=l ,7=1 

The coupling constants are normalized in a convenient way so that the derivatives are 

dl+l 
K=l 

We may define two sequences of mutually orthogonal monic polynomials 7r„ (x) , cr„ (y) of degree n such 
that 

/ da; / dy7r„(a;)a™(y)e-^i(^)-^^(^)+^^ = /i„5„„ , (2-3) 

TTn{x) ^x'" + ... , a„(y) = y" + ... . (2-4) 

In order that the integrals be convergent, one should suitably define the two closed contours of integration 
Tx,Ty. If we require that these be the real axis, the degrees of the potentials must be even with the 



8 



leading coefficient having positive real part. In applications of random matrices to string theory however, 
the integral is not convergent on the real axis, and the contour should approach oo in some appropriate 
Stokes sector in the complex plane. 

It is more convenient to deal with the quasi-polynomials defined by 



1 



/En 



Tr„{x)e' 



-Vi{x) 



and their Fourier-Laplace transforms 



1 



o'n(2/)e' 



-V2iy) 



dye-yMy) 



(2-5) 



(2-6) 



The choice of normalization is somewhat arbitrary. We have here chosen the most "symmetric" one, 
in which the leading coefhcients in the various recursion relations are the same. In this notation the 
orthogonality relations take on a simpler form. 



da:V'n(a;)^^(a;) = J dy'ip^^{y)(l),n{y) ^ J J dxdy>pnix)(j)miy)e''y 



(2-7) 



(We suppress for the present the specification of the contour of integration.) We shall think of the spaces 
spanned by the ?A„(a;)'s and (/)„(2/)'s as infinite graded spaces in duality through the pairing in eq. (2-7). 
It can easily be seen, using these relations and integration by parts, that multiplication of the ipnS by 
X produces a linear combination of tAm's with n — d2<m<n+l and multiplication of the 0„'s by y 
produces a linear combination of 4>j 's with n — di<j<n + l. Moreover it is clear (through integration 
by parts) that multiplication of the ^„'s by x is dual to application of dy to the (/)„'s and vice- versa. 



2.2 Recursion relations and generalized Christoffel— Darboux formulae 

We denote by Q and P the semi-infinite matrices which implement multiplication and differentiation by 
X on the space spanned by the ipn quasi-polynomials. Introducing the semi-infinite column vectors 



* [-00, ■ 

oo 

the above remarks imply that 



and $ 



X \I/ 

oo 

y $ 



oc 



d 

— * := 

dx oo 

d 

— $ 

dy oo 



where P and Q are semi-infinite matrices of the form 



Q:= 



ao(0) 
«i(l) 



7(0) 
ao(l) 





7(1) 



OO 



o:d2{d2) ad2-i{d2) 




"0(^2) 7(^2) 



(2-8) 



(2-9) 



(2-10) 
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-p 



satisfying tlic string equation. 



/3o(0) /3i(l) ••• PdAdi) 

7(0) /3o(l) /3i(2) ■•• + 

7(1) /3o(2) 

7(2) /3o(3) ■•• 



(2-11) 



(2-12) 



The fact that both matrices have a finite band size as indicated in (2-10), (2-11) follows from the fact 
that they are related polynomially to each other through the potentials Vi and V2 as follows. 



Lemma 2.1 The two matrices P and Q satisfy the following relations 

(p+v;(o))^„ = o , 

(-Q* + V2'(-P*))^„ = , 

where the subscript >o means the part above the main diagonal (main diagonal included). 
Proof. It is obvious that 



(2-13) 
(2-14) 



(F + y/(Q))vI/ 



hi+1 



^Jn+i{x) + lower terms. 



-^ + Vi{x)]M^) 



<(a;)e" 



-Vi{x) 



-tpn-iix) + lower terms. 



Eq. (2-15) means that Q has only one diagonal above the main diagonal with entries given by 



7(n) = \/hn+l/hn 



(2-15) 



(2-16) 



(2-17) 



Eq. ( ^-16| ) then implies that eq. ( 2-13| ) holds and P has di diagonals above the main one. Repeating 
the argument for the 0„ quasi-polynomials similarly shows that — P* (i.e. multiplication by y) has one 
diagonal above the main one, eq. (2-14) holds and — Q* (i.e. differentiation by y) has d2 upper diagonals. 
This proves the lemma and also shows that P and Q are of the finite band sizes indicated in ( ^-10| ), 
(p^). Q.E.D. 



Eqs. (2-£) are just equivalent to the following set of recursion relations 

d2 



00 

P^ 



xtpn{x) = 'y{n)ipn+i{x) + ^aj{n)'tpn-]{x) 

d '^^ 
— ■0„(a;) = -lin- l)V'„_i(x) (3 ;j(n + 3)^:^+0 {x) 



(2-18) 
(2-19) 
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di 

y4>n{y) ^ 7(n)(/)„+i(y) + 

>n(y) = -7(n-l)(/'„-i(y)-^aj(n + j)(/>„+j(y) 
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(2-20) 
(2-21) 



In the following we also define 

a_i (71) :=7(n) (n) , a,(n) 0, ^ [-1, ^2] , /3,(n) := , Vj ^ [-1, rfi]. (2-22) 

Defining similarly the semi-infinite row vectors consisting of the Fourier-Laplace transformed functions 



(2-23) 



it follows from the dual pairing (2-7), and integration by parts that: 
Lemma 2.2 



00 
00 

Now, we introduce the shift matrices 

10 

10 

A:= 1 

1 





i=-i 

d2 



i = -l 



(2-24) 
(2-25) 
(2-26) 
(2-27) 



i=-i 



:= A* 





1 
10 
10 

■■• 



(2-28) 



(The notation A^^ is a convenient shorthand for the transpose A*, but only signifies that A~^ is the 
right inverse of A. It leads to the abbreviated form A^-' := (A*)-'.) Introducing the diagonal semi-infinite 
matrices 

a, diag(a,(0),a,(l), . . .) , /?, := diag(/3,(0),/3,(l), . . .) , (2-29) 
(where we set aj{n) := and /3j(n) :— when n < j), eqs. ( 2-10 ) , ( ^-11 ) can be concisely written as 



d2 dl 



(2-30) 
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The commutation relation in eq. (2-12) gives in particular the following quadratic relations between 
the coefficients {aj, (3k}, 

di 

(2-31) 



di 
3=1 


di 

= ^ {n)Pj (n + l + l) , 

3=1 


yn,yie [0,di] 


d2 

Pk~i~-i{n + k)ak {n -f k) 

k=l 


d2 

= y^^akjn + I + l)(3k-i-i{n) , 

k=l 


Vn, V/ e [0,d2 



(2-32) 

Our first objective is to define a set of closed differential and difference systems for vectors con- 
sisting of finite sequences of the functions {ipn}, {0n} and their Fourier-Laplace transforms, equivalent 
to the systems ( 2-l^ )-(2-21), and to study their properties and relations. Consider first the sequence 
of functions {'(/'„(x)}. For these we have a multiplicative recursion relation defined by the coefficients 
{aj{n)} and a differential recursion relation defined by the coefficients {(3j{n)}. We shall show presently 
that we can define closed systems of first order linear ODE's for any consecutive sequence of ^2 + 1 
functions {ij}N-d2i ■ ■ ■ i V'iv) (or ('/'jy-i' ' ' ' ' ^N+d -0^ with coefficients that are polynomials in x. Similar 
systems can be constructed for any sequence {<pN-di, • ■ • , (('n) (or (V'^.^i ■ • ■ ; '^^ N+d introduce 
the following definitions and notations 

Definition 2.1 A window of size di or d2 is any subset of di or d2 consecutive elements of type tpn, 
4> , (pn or ip , with the notations 



N 



* := [i/.iv,...,V^ArH^i]* , N>0, $ := [(/.Ar,...,0jvw.]* , A^>0 (2-34) 



N 
N 



N>d2, 


$ 




N 




N 


TV > 0, 


$ : 


N>d2, 


$ 




N 




N 


N>0, 





I* , N>di (2-33) 



*^=[^^^,'--->J ' N>d^^ N>di (2-35) 



^■■=itN^---^tN^,]^ N>0, $:=[0^,...,^^^J, A^>0. (2-36) 

Notice the difference in the positioning of the windows for the vectors constructed from the ipnS and the 
4>nS, and the fact that the barred quantities are defined to be row vectors while the unbarred ones are 
column vectors. 

Definition 2.2 For any N for which these are defined, the pairs of windows (^vj/, ^ ) as well as ^ ) 

N N 

of dimensions 1^2 -f 1 and di -t- 1, respectively, will be called dual windows. 

The reason for identifying these particular windows as dual will appear in the sequel. 
Let us now consider the kernels 

iV . ^1 . N 



Kii{x,x') := Mx)^Jx') , K22{y',y) ^ ±Jy')My) , (2-37) 

71=0 n=0 

Ki2{x,y) }2 Mx)4>n{y) , K2i{y\x') ^ ±Sy')tS^') ' (2-38) 



n=0 n=0 

that appear in the computation of correlation functions for 2-matrix models |22[ 
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Define the following pair of matrices, which will play an important role in what follows: 



N 



N 



U 


u 




U 


U 










a2iN) 


aim 










-1) 




ai{N + l) 































ad,{N+d2-l) 



















^(iV-l)1 


f3dM 








/3i(A0 














Pi{N + l) 













Pd,iN+2) 
























(2-39) 



(2-40) 



For any N , the r ecursion relations (|2-18|) , ( ^-20| ), (|2-24|) , ( ^-26| ) and the differential relations ( ^-19|) , 



(2-21), ( 2-25 ), ( 2-27 ) imply that the following generalized Christoffel-Darboux formulae, as well as their 
"differential" analogs are satisfied. 

Proposition 2.1 Generalized Christoffel-Darboux relations: 



N 



{x - x)Kii{x,x) = 



N 



{y' - y)K22{y' ,y) = 



d2 J-1 

7(iV - l)^Ar0^_, - E E + >')^N+k'J'N+''-l 

j=l k=0 

'N~l N \ 

$ (x'),A^'(a;) , 

N J 

di j-l 

N-j+k 

j=l k=0 

N-1 N ^ 

^ (y'),B$(2/) 

N 



(2-41) 



"Differential" generalized Christoffel-Darboux relations: 

{d^' +d^)Kii{x',x) = - f $(x'), (B)* V(x) 



N 



{dy, + dy) K22{v\ y) - - My'). {aY $ (y) 



N , N-1 



(2-42) 

(2-43) 
(2-44) 



Proof. Use the relations ( 2-18| )-( ^-21 ), ( 2-24 )-( 2-27 ) and simplify the telescopic sums by cancellation of 
common terms. Q.E.D. 

Although it will not be needed in the remainder of this paper, for the sake of completeness, we also 
include the following analogous result for the kernels K12 and K21, which may be similarly derived. It is 
related to the above by applying Fourier-Laplace transforms with respect to one of the variables. 



Proposition 2.2 



N 



'N-1 



ix + dy)Ki2ix,y) 



N 



iy),A^ix) 

N 



(2-45) 
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N 



'N-1 



N 



N 



[y' -d.,,)K2i{y',x') - 



N 



{dy,-x')K2i{y'.x') 



(y'),B|(2/)j 

N-l N \ 

B$(a;') , 

N J 

N-l N 

$ ix'),A:^\y' 

N 



(2-46) 

(2-47) 
(2-48) 



2.3 Folding 

We now intr 
respectively. 



We now introduce the sequence of companion-like matrices a(x) and h{y) of sizes d2 + I and di + 1, 

N N 



N 








7(7V) 


















1 

-ai{N) (x-ao(jV)) 
fiN) j{N) 




1 

-l3i{N) (y-l3o(N)) 
-riN) jiN) 



N>d2 



(2-49) 



N>di 



(2-50) 



We then have the following: 



Lemma 2.3 The sequence of matrices a, b implement the shift N ^ N + 1 in the windows of quasi-polynomials 

N 

in the sense that 



and in general 



a vl>(x) = vj/ (x) 

N N N+1 



xj/ = a • ■ a xj/ , 

JV4^- JV+i-1 iV iV 



b$(y) = $ (y) 

N N N+l 



$ = b ■■■b$ 

N+j N+j-1 NN 



(2-51) 



(2-52) 



Proof. This is nothing but a matricial form of the sequence of recursion relations ( 2-18 ), ( 2-20 ) expressing 
the higher order polynomials as linear combinations of a fixed subset with polynomial coefficients. Q.E.D. 

We will refer to this process of expressing any ipn (x) by means of linear combinations of elements in 
a specific window with polynomial coefficients as folding onto the specified window. 



The determinants of the matrices a and b are easily computed to be 

N N 



det(a) = (-l)'^^+V,,(7V)/7(^) 



From eqs. ( [2-13| ), ( 2-14 ) we find the relations 



det(b) = (-l)'^^+V.,(^)/7(A^) 

N 



(2-53) 



(2-54) 



Since the coefficients 7(A^) are the square roots of the ratios of normalization factors, they cannot vanish 
for any N, and neither can ad2{N) or Pdi{N), since the deformation parameters Ud^+i, Vd2+i are the 
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leading coefficients of the polynomials Vi{x), V2{y) and hence also may not vanish. It follows that the 
matrices a and b are all invertible. We denote their inverses as follows 

N N 



N 

a 



N 

b 



-aa^i{N) 
1 




-)3di-i(jV) 
/3<ii(JV) 

1 










x-ao{N) -7(jV) 







1 

V-0o(n) 












(2-55) 



(2-56) 



The shifts ^ — 1 arc thus implcnK^ntcd by the inverse matrices a, b, and the folding may take place 
in either direction with respect to polynomial degrees. 



2.4 Folded lineeir differential systems 

We now define the following sequences of finite diagonal matrices 



N 



aj:=di&g[aj(N+j-di),aj{N + j-di + l),...,aj(N+j)] , j = -1, . . .d2 , (2-57) 



N 



Pf= diag [/3, {N + j- d2), Pj{N + j-d2 + l),..., 13 j {N + j)] , j = -l,... di 



(2-58) 



N 



Recall that a_i(n) = 7(n) = /3_i(n) by our previous conventions, but the diagonal matrices a_i and 

JV 

differ in dimensions. When the context leaves no doubt as to the dimension we will write them as 



or 



a_i=7:= diag [^{N - rfi - 1), . . . , 7(iV - 1)] 



2_i=7:= diag [7(iV - d2 - 1), l{N - da), . . • , 7(^ " 1)] • 



In either case, we denote the inverse matrix as 



N 



7 := (7)"^ . 

JV 

We can now give the closed differential systems referred to previously. 

Lemma 2.4 The windows of quasi-polynomials ^, $ satisfy the following differential systems 

N N 



(2-59) 
(2-60) 

(2-61) 



— * 

dx N 



-Di{x)-iS/ , N>d2 + 1, 

N 



Q N 

— $ = -D2{y)^ , A/'>di + l, 
ay N N 



(2-62) 
(2-63) 
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where 



N N N-1 ^ ■sA ^ 

Di{x)--1 a + Po + } Pj a a ■ 

^ ' ^ N+j-lN+j-2 

J = l 



AT 



d2 



±- , , N N-1 iv \ - 

D2{y) :=7 b + ao +}_^ 



N 



N+j-lN+j-2 



■a G gld2+ilx] 



■b e gldi+il 

N 



(2-64) 
(2-65) 



These , taken together for all N are equivalent to the relations (2-li), {2-21) when the recursion relations (2-51), 
(2-5i) are taken into account. 



Proof. Consider the case of the ■i/'n's- The differential relations (2-19) may be written by stacking them 
in a window of size (i2 + 1 , as follows 



d 



dx N 



N 

AT AT 1 ' ^ AT I ^■ 



JV-1 



N+j 



(2-66) 



Using the folding relations e g. (p-52 ), we immediately obtain (2-62), ( 2-641 ). The same procedure applied 
to eq. ( p^ ) yields ( |2^ , (|2-65D. Q.E.D. 
We can repeat a similar procedure for the sequences {ip iy)}neN ^^id {0 (x)}„£n. The corresponding 

AT N 

windows are represented as row vectors ^ and $ since their components are naturally dual to the (j)n's 
and ipnS respectively. The matrices defining the relevant folding are now 



N 

a 



N 

b 



7(Ar-l) 

-ai(Af-H) 
7(JV~1) 



-ad3(jV+rf2) 
7(JV-1) 

7(Af-l) 



-/3rf^(jV+di) 
7(Ar-l) 



1 











1 



1 
'•• 
1 




(2-67) 



(2-68) 



and we again denote their inverses as 



a 

N 



b := 

N 



(2-69) 



As previously , we now have: 

.AT iV 

Lemma 2.5 The sequence of matrices {a} {b} implement the shift N i 



iV ■ 



JV-l AT AT 

* =*b 



AT-l AT ^ 

$ — $ a. 



(2-70) 
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N 



N 



N 



N 



Similarly to the diagonal matrices aj , (3j , we define the matrices a„ and (3 . as 



N 



a, :^ dia,g{aAN),aj{N + 1),. . .,aj{N + di)) , j 



:= diag (N), /?, (iV + 1), . . . (3, (N + d^)) , j 



(2-71) 
(2-72) 



As before we have the two definitions 



a_i 7 := diag (7(iV), j{N +l),^...j(N + d,)) 



N 



N 



■■= 7 diag (jiN), jiN +l),...,j{N + d^)) 



(2-73) 
(2-74) 



which will be used if there is no ambiguity regarding dimensions. 

By repeating a procedure similar to what led to the differential systems in Lemma 2.4, we find: 

JV-l JV-l 

Lemma 2.6 The dual windows of Laplace-transformed quasi-polynomials ^ , $ satisfy the following differ- 
ential systems 



where 



Q N-1 

d-y^^y^ 



N-1 N 

* {y)D,,{y) , N>di + 1, 



r) N-1 N-1 N 

— $(x) =^ix)D,{x), N>d2 + 1, 



(2-75) 
(2-76) 



11 iV-l N-1 r-^iv-j.jv- 

My) :=b7+Qo+2^b b 



JV-l N-2 N-j pf_i 

b 



AT 

^l(^) 



dl 



N-1 - ■r-^JV-lJV-2 iV-j-'V-l 



Summarizing, we have thus obtained four differential systems 



Size (d2 + l)x {d2 + 1) 



d ^ 

— m{x) = -Di{x)-^{x) 

ox N N 
g N-1 N-1 N 

— $ (x) = $ {x)D,{x) 



Size (dl + 1) X {dl + 1) 



g N-1 N-1 N 

oy 



9 , . 
ay N 



N 

D2{y)Hy) 

N 



(2-77) 
(2-78) 



(2-79) 



It should be noted that the two matrices Di and (as well as D2 and D.2) have so far only superficial 
similarities. In particular they do not depend on the same subsets of the coefficients {aj{n)} and {f3j{n)}. 
On the other hand the pairs {Di^D_2) and (I?2,i2i) depend on the same aj(n)'s and f3j{n)^s although 
they are of different dimensions. 



2.5 Deformation equations 

The following Lemma gives the effect of an infinitesimal deformation in the coefficients {uk, vj} expressed 
as differential equations for the (semi)-infinite vectors xpooj $00 of biorthogonal quasi-polynomials (as well 
as their Fourier-Laplace transforms) and for the matrices P, Q. (Derivations with different conventions 
can be found in |l|, |, |l|.) 
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Lemma 2.7 





* = 

oo 


oo 






* = 

oo 




oo 




oo 




)'* 

oo 




$ = 







where 



K 

Componentwise these read, 



where we have used the notation 



J 
if 

J 

dvj(l}u{y) = ^ V'/(n)(j!>„+j(y) , 



C/f (n) := (7^,.+„ , T//(n) := 



(2-80) 

(2-81) 

(2-82) 
(2-83) 

(2-84) 

(2-85) 
(2-86) 
(2-87) 
(2-88) 

(2-89) 



(The same relations hold for the Fourier-Laplace transforms with respect to the variables x and y, since the 
coefficients do not depend on these variables.) 

Moreover, we have the following equations for the matrices P, Q: 



(2-90) 
(2-91) 

du^P = -[P,U''] , (2-92) 
d^jP ^[P,V''] . (2-93) 

Proof: Equations ( ^-80D and ( |2-83D are just definitions, ( |2-82| ) and ( |2-8lD follow from ( |2^ . Eq. ( |2-84D 
is proved in a way similar to Lemma 2.1. From the definitions (2-5), one has: 



and 



2 K 



'hn 



K 



Since 9u^7rn(x) has a degree lower than n, one sees that eq. ( 2-94 ) impUes that 



(2-94) 
(2-95) 

(2-96) 
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and similarly eq. ( ^-95| ) implies that 
They also imply that the diagonal part must be 



1 hn 1_ t(^K\ 

2 hn ^ 2K^^ ' 



(2-97) 



(2-98) 



which proves ( 2-! 

The equations (^-90|)-(p^) follow from multiplying ( |2-80D and ( ^-8l| ) by x and ( |2-82| ) and ( ^-83| ) by 
y and the linear independence of the component functions forming the vectors The coefficients of the 
expansion must vanish, and that is precisely the relations (2-90) to (2-93). Q.E.D. 

We also require the "folded" version of the deformation equations. This leads to eight equations giving 
the action of dvj on vj/, J[. , $ and $ . They are introduced in the following lemma, for which 



N 



N 



we ne e d to defi ne di agonal matrices which play roles similar to that of the matrices defined in ( ^-57 ), 
( ^-58| ) ( |2-7lD , ( ^-72D for the differential equations with respect to a; or y. 



N,d 



U j,K cliag(C// (7V-d), . . . ,Uf {N)) , 
V ,,j diag(l^/(7V-d-j), . . . , y/(iV-i)) . 

With this notation we have: 

Lemma 2.8 The deformation equations can he written in the folded windows (and dual windows) as 

g jv-i jv-i AT,* 



J. *= Uk* , 

OUk n n 

d 

— * = - V J* 

OVj N N 



7^ * = * , 
OUk 

P) N-l N-1 iV.* 

— ^ = - ^ V J , 

OVj 



(2-99) 
(2-100) 



(2-101) 



where 



Q N,'S 

t: <& = - U K$ 

OUk n n 
Q N,<S> 

= V j<E> , 

OVj N N 



V,* . , N,d2 

TJ K ■- } U j,K a ■ ■ ■ a 

EN,d2 N-j N-l 
V j,j a • ■ ■ a , 

JI, N-j,di N-j N-l 
E U j.K h ■■■ h 

j=0 

N+j.di 



JV,* 

V J 



V J.J b ■■■b 



g N-l iV-lJV,« 
$ = - $ u 



a — ^ ^ if I 

OUK 
f) N-l N-1N,<S> 

7^ t = * , 
OVj 



JV,* 

U , 



iV.* 

Y , 



Uk 



N,<S> 

V , 



K 






• b 




N+j-l 


N-l 

E b ■ 


N-J 

h 






K 




E ^ 


■ ■ a 






J 




Es- 


a 




N-l+j 



N+di-l,di 

U j,K , 



V 



u 



N+j+d2-l,d2 

V j,J 



(2-102) 



(2-103) 



(2-104) 



Proof. This follows exactly the same lines as the proof of Lemma 2.4 



Q.E.D. 
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3 Compatibility of the finite difference— differential— deformation 
systems 

We want to now prove that the recursion relations (2-51), ( p-7C ) , the hnear differential systems (2-79) 
and the systems of deformation equations (2-101), ( p-102| ) are all compatible in the sense that they admit 
a basis of simultaneous solutions (fundamental systems). 

Proposition 3.1 The shifts i-^ + 1 in eqs. ( ^-5l| ) implemented by a and b and the sequence of 
differential equations ( P-62| ), ( |2-63D , respectively, are compatible as vector differential-difference systems. 
That is, there exists a sequence of {d2 + 1) x (c?2 + 1) fundamental matrix solutions < 'if{x) > and 



N 



N>d2 + 1 



{di + 1) X (c?i + 1) fundamental matrix solutions < $(y) > simultaneously satisfying the equations 



N 



N>di + 1 



^{x)^a{x)^{x) , (3-1) 

N+l N N 

d N 

^^{x) = -Di{x)^{x) , N>d2 + l, (3-2) 



dx N N 

and 

* (y)-b(?y)*(y) , (3-3) 

N+l N N 

d N 

— ^{y)^-D2{y)^{y). N>di + l, (3-4) 
oy N N 

respectively. The same result holds for the barred quantities and the shifts N ^ N — 1 implemented by 
w w . . . 

a and b. That is, there exist fundamental solutions { ^ {y)}N>d2i+i of dimension {d2 + 1) x [d2 + 1) 

N-l 

and fundamental solutions { J& (a;)}Ar>d2+i of dimension [di + 1) x {di + 1) simultaneously satisfying the 
recursion relations and differential systems 

N-l N N 

* (y) = *(?/) b(y) , (3-5) 



dy 

and 



f) N-l N~l N 

TT * (y) = * iy)D2{y) , iv > rfi + 1 , (3-6) 



Af-1 N 



* (x) = *(x) a(a;) , (3-7) 
d_ 

dx 



f) N-l N-l N 

^ {x) ^ ^ {x)D^{x) , N>d2 + 1, (3-8) 



respectively. 



Proof. We prove the compatibility for only one of the four shift-differential systems, the others being 
completely analogous. 

The statement amounts to proving that 

N TV / W+1 \ N+l N d 

d^+Di(x) = ao + Di(x) o a = + a(a:) DAx) a(j;) + a(x)— a(x) , (3-9) 

\ I N N dx N 
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where the dependence of a on a; has been emphasized. 

N 

Let 



^{x):=\ijN^,{x),...,iPN{x) , N>d2 (3-10) 



be any solution to the equation 



9^ + Di(a;) ) |(x) = . (3-11) 



At this stage the labeUng N — d2, . ■ . ,N has no particular meaning because there are no Vn(a;)'s with 
n ^ [N — d2, N]. Nevertheless we can define 

^ (x) := a • • • a *(x) j>l , (3-12) 

and 

^ {x):=\^ ■■■^a^^ix) , 0<j<N-d2 (3-13) 

JV-^ N 

Because of the recursive structure of the matrices a, the above definition, e.g., for ^ actually defines 

not d2 + 1 new functions, but only one new function: V'AM-i- Therefore, componentwise, we have defined 
a sequence of new functions ipm (x) , which satisfy the recursion relation 

X'ipm{x) = ^ aj{m)'ipm-j{x) ,m > ^2 • (3-14) 

By this definition and by the structure of the matrix 

N N N-1 N ^ N 

Di =7 a -h /3o + > /3j a a • • • a , (3-15) 

the differential system componentwise now reads 

di 

d^Mx) = - ^ /3j(n + j)i>n+j{x) , n = N-d2,...,N , (3-16) 

where the ^„'s that fall outside the window N — d2, ■ . ■ ,N have been defined above in terms of the ones 
within the window. Therefore we need to prove that the newly defined function 

-• x-ao{N)~ J^ai{N)~ 

satisfies the same sort of diff'erential equation as the preccdiug ones. (A simple argument by induction 
then shows that all ^pN+j satisfy the same sort of difi^erential equation for any j > 1). This in turn 
amounts to proving that 

di 

d^xl>N+i{x) = - ^ /3,(iV + 1 N+i+j{x) . (3-18) 
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To do this we compute 

j{N)dj:'lijN+l{x) = 



dx 



xiPn{x) - ^ ai{N)tpN-i{x) 



1=0 



d2 di di 

2 J2 ai{N)pj{N-l+j)^N-i+o{x)+^N{x)~x J2 PjiN + MN+j{x) 
1=0 j=-i j=-i 

d2 di 

Oil{N)f},{N~l+j)i^N-l+o{x)+i>N{x) 

1=0 j=-l 
d2 di 

-JlJl ai{N + j)P,(N + j)i:N+,-i{x) . 
i=-i]=~i 



(3-19) 



Rearranging eqs. (3-18, 3-19), we have to prove the identity 



i=-i 



d2 di 

Y J2 aiWPjiN^ l+j)^N-i+j{x) + i^N{x) 

l=Oj=-l 
d2 di 

- E E + J)^^- + j)^N+j-i ix) , (3-20) 

i=-ij=-i 



or equivalently 



d2 di 



i^N^x) 



- J2 E aiiN)Pj{N~-l+j)^M^i+,{x) 
i=-ij=-i 

d2 di 

- J2 E MN + .miN + j)iJN+,^iix) . 
i=-ij=-i 



(3-21) 



But this last equation is nothing but the Heisenberg commutation relations in eq. (2-31) and eq. ( ^-32| ) . 
This means that rearranging the coefficients in front of il)M+r(,x) in the RHS of eq. (3-21), the only 
nonvanishing coefficient is that of ?/'Ar(x) and it is exactly 1. A similar argument may be used to prove 
that the relations (3-16) hold also for 1 <n < N — d2 

For future convenience we remark that this verification amounts to the fact that the coefficients of 
the ipN+rixYs are the same as for the orthogonal quasi-polynomials ij^N+rix), since it relies only on 
the recursion relations, which are the same. Given that the quasi-polynomials TpN+rix) are linearly 
independent, the equality of LHS and RHS follows also for any other sequence. Q.E.D. 

Prop. 3.1 means that we can define d2 + I sequences of functions | xjjn '^ (x) \ in such a way 

L J neN,q=0...d2 

that in any "window" of size d2 + 1 they constitute a fundamental system of solutions to the differential 
system (3^). One of these sequences is obviously provided by the orthogonal quasi-polynomials. Each 
of them satisfies both the recursion relations 



i=-i 



and the derivative relations 



d^^Pi^Hx) = - J2 f3,in+Ml:l,{x) 



(3-22) 



(3-23) 
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Remark 3.1 In principle, in order to define these d2 + 1 sequences one should solve the differential system ( B-2) 
in a given window and then define recursively the rest of the sequence backwards and forwards. To pass from 
the semi-infinite case to the infinite one, we may define the full sequence t/)!'' for n € Z just by application of 
products of the matrices a and their inverses, provided the aj(n)'s are so defined that all the a's are invertible. 



In a completely parallel manner we can define di + 1 sequences of functions I (pn ^^ (y) \ which 

L J neN,q=0...di 

provide fundamental systems satisfying 



dy + D2iy) 



My) = . (3-24) 

N 



Moreover, with minor modifications, we can construct analogous sequences for the dual systems 

N-l N-1 N 

5,* (2/) = *(2/)^2(y) (3-25) 

N-l N-l N 

9, *(x) = ^{x)D,{x). (3-26) 

The only difference is that the matrices a and b now implement the shift iV ^ iV — 1. The barred 
sequences I 6^'^^ (y) > and I tb^'^^ (y) > will therefore satisfy the recursion relations 

- E"'(" + 0£t(^)' ^-^n^(^) ^ f3An)£%ix) , (3-27) 

l=-l J=-l 

vi^niy) - E /5.(" + j)£t (^) ' ^y'l^n^y) = E «K"Mi'^l,(y) ■ (3-28) 

j=-l l=-l 

A completely analogous statement holds for the matrices defining the deformation equations in any 
window. 



Proposition 3.2 The shifts N N + 1 implemented by a in eq.(3-l) and the sequence of differential 

N ^ ^ 

equations 

N,-9 AT,* 

5„^*= Ux*, 5„,* = -Vj*, N>d2 + 1 (3-29) 

N N ' N N 

N 

are compatible, as are the shifts N N — 1 implemented by b in (|3-5| ) and the sequence of differential 
equations 

* = * U ^, 9.., * = - * Y , N>1 . (3-30) 

N 

Similarly the shifts implemented by b and a are compatible with the equations 

N 

Af.$ AT,* 

* = - U K * , (9^, * = V J * , N>di + l (3-31) 

N N ' N N 

N-l N-1N,<S> N-l N-1N,<S> 

du^- ^ = - * U ^ , * = * Y J , N>1 . (3-32) 
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Proof. We will prove compatibility of only one of the eight kinds of systems with the shift; the remaining 
cases are proven similarly. 

As for Prop. 3.1, we first define a continuous parametric family (depending on x) of solutions to the 
system 

TV,* 
du^- - U K 



^ = Q 

N 



(3-33) 



We then define the shifted functions := a • • • a vl/ so that the equation reads, componentwise 

N+j N+j-1 N N 



K 



duK'4'n = U^{n)'ij}n+j , 

3=0 



= N -d-2 



(3-34) 



Then we have to check that the newly defined ipN+i also satisfies 

K 

du^^N+i = J2 Uj''{N+l)i'N+i+j (3-35) 
j=o 

(and by induction the corresponding equations for ^^^r, f As before, we use the relations 



X — ao{N) ~ 



d2 



ai{N) - 



IpN-lix) 



(3-36) 



and the differential system satisfied by the 'ipN~d2j ■ ■ ■ i V'JV- To conclude the equality we can reason as in 



the remark in the proof of Prop. 3.f , Q.E.D. 

The final proposition in this section assures that the deformation equations are compatible with the 
x,y differential equations in all windows. 

Proposition 3.3 The system of equations 



N 

N 



N.l/ 

duK - V k] *(2^) 



TV.* 

9vj + V (J) 



N 
N 

* (x) 

N+1 



ai{x)^{X), N>d2 + 1 

N N 



is compatible, and hence sequences of fundamental systems of solutions < '^(x) 
exist. The same statement holds for the system 

N 



N>d2 



dy + D2] My) 



N,'S> \ 

J N 

N.<!> \ 

9vj - V (J) j My) 



* (y) 

N+l 



, 
, 
, 



h{x)^{y) , N>di + 1 

N N 



(3-37) 

(3-38) 

(3-39) 
(3-40) 

to all equations 

(3-41) 

(3-42) 

(3-43) 
(3-44) 



24 



The corresponding systems for the barred sequences are also compatble, and hence also admit simulta- 

N-l N-l 

neous sequences of fundamental solutions ^ (x) and J& (y). 

Proof. The proof will only be given for the system (|3-37|) -(F40|) since the others are proved in the same 
way. The comp a tibilit y follows from Props. and 3^ together with a proof of compatibility of the 
equations ( 3-37 )-( 3^^ ). Indeed, from the ^2 + 1 functions in 



* = [ijN-d2{x,U,v), . 



, ,-0Ar(a;,u,w)]* 



(3-45) 



we can consistently define a whole sequence of functions V'n's by means of the x— recursion relations in 
such a way that componentwise the system reads 



di 

E 



n+j , 



K 

J 

3=0 



Taking cross derivatives and using these expressions one gets, e.g., 

K di di 



f 9 



k=Oj=-l 
di K 

■ E Y.^k{n)l3j{n + j + k)ijr,+j+k 

j=-l k=0 



\duK 



Pjin + j) ^. 



(3-46) 
(3-47) 
(3-48) 

(3-49) 
(3-50) 



The e xpressions for the derivatives of the coefficie nts (3 j may be obtained from the deformation equation 
( 2-92 ) for P and then substituted back into eq. ( 3-49 ). However, to prove the equality of the two cross 
derivatives it is sufficient to collect the coefficients of i^n+q and note that exactly the same coefficients 
appear when the functions {ipn} are replaced by the orthogonal quasi-polynomials {ipn\, for which the 
equality of the two expressions certainly holds. Since the orthogonal quasi-polynomials are linearly 
independent functions, the individual coeficients must agree. (This is essentially the same argument as in 
the remark at the end of proof of Prop. 3T). The mutual compatibility of the du^ and d^ , d eform ations 
is proved in exactly the same way. One just takes the the cross derivatives in equations (3-47) and (3-' 



and notes that, since these are the same as for the case of the orthogonal quasi-polynomials {ipn\, the 
corresponding coefficients in the cross differentiated expression must be equal. Q.E.D. 



4 Spectral Duality 



The aim of this section is to st ate and prove some remarkable relations between systems related to dual 
windows (in the sense of Def. 2T), which will justify the terminology. One of the main results will be 



that the four spectral curves given by the characteristic polynomials of Di,D^, D2,D.2 associated to the 

N-l N-l 

four systems (eq. ( 2-79 )) on the two pairs of dual windows (\I>, $ ), ($, ^ ) are actually the same curve. 



N 



N 
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4.1 Dual spectral curves 

First we need a linear algebra lemma. 

Lemma 4.1 Let T be a square matrix having the block form 






Fi 

















F2 


































Fd 




Gi 


Ga 




Gd 



where the d + 1 blocks have compatible sizes and the diagonal blocks are square. Then 

det [1 - T] = det [1-D] , where 

d-l 

D := Gd + ^2 • Fh+i ■ ■ ■ Fd , 

fc=0 

and 1 denotes, according to the context, the unit matrix of appropriate size. 
Proof. Let 



F ■- 



' 


Fi 








" 








F2 


































Fd 


. 















We multiply the matrix 1 — T from the right by the matrix 



1 


Fi 


FiF2 


FiFiF.i 


Fi 


■■■Fd' 





1 


F2 


F2F3 


F2 


■■■Fd 























1 


Fd 














1 



(4-1) 



(4-2) 
(4-3) 



(4-4) 



(4-5) 



Since the matrix (1 — Np) ^ is unimodular, the determinant of 1 — T remains unaffected. Then one 
computes 



{1-T)-{1-Nf)-' 



1 











" 





1 

















1 






























1-D 



from which the statement follows by taking the determinant. 



(4-6) 



Q.E.D. 
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N N N N 

Proposition 4.1 The spectral curves associated to the characteristic polynomials of Di, D,2, D2, D.i 
are pairwise equal. More precisely, we have the formulae 



Udi+i det 
Udi+i det 



xl 



N 



di + 1 



D2{y) 



N 



= Vd^i det 
= Vd^i dot 



N 



yld2+i -~D-i_{x) 



N 



yld2+i - Di{x) 



(4-7) 
(4-8) 



which connect the spectral curves of the differential operators of different dimensions operating on the 
two pairs of dual windows. 

Proof. We will only prove one equality since the other is proved similarly. We start with the computation 
of the characteristic polynomial of Di [x) 



det 



yl-D^{x) 



det 



N N-1 

7 a 



det 



/ N ^ N 
1 - a 7 yl- 00- y 0i a • • • a 



Vd2 + 1 



det 



i=i 

N 



/ N ^ N 



a •a 

N+j-l N 



where we have used the identity in eq. (2-54). Now we use Lemma 4.1 with the identifications 



N 

Gd,= a 7 ( yl- 00 



N 



and thus obtain 



det 



N 



N 



1 - a 7 yl- 0o-yZ0j a • • • a 

N-l N \ ^ N+j-1 N 

The matrix Tab is defined by 



det 



+ I)(d2 + 1) - Tab 



T, 



ab 






a 




























a 

N 


N 

-a /3di7 


N 

~SL /3di_i7 

N-l N 




N 



a 








" 


N+di-l 

























a 









N 













a 








Af-l. 






1 




























1 


N 

-0d^l 

N 


N 

-/?di-i7 

N 




N 

(2/I- /3o)7 



(4-9) 



(4-10) 



(4-11) 



(4-12) 



(4-13) 
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We regard Tab endomorphism of C^^^^ (E) C'l^"'^. Let P12 be the involution interchanging the two 

factors of the tensor product and let the matrix C implement the reversal of order endomorphism within 
the (^2 + 1) X {d2 + 1) blocks: 



d2+i times 
C Blockdiag {R, R,..,R) , R := 



" 











1 " 











1 






















1 











1 















(4-14) 



A direct inspection shows 






b 


































JV-2 

b 

















N-l 

b 


N N-l 
-b "da 7 

N-l 


N N-l 

-barf,_i7 

Af-l 




Nn-1 

-btti 7 

N-l 


b (xl- ag J 7 

N-l 



(4-15) 



where we now have (i2 + 1 square blocks of dimension di -I- 1 (i.e. the number of blocks and the dimensions 
of the blocks have been interchanged). The barred symbols are precisely those defined in eqs. (|^,|j|), 

(FtiI, Fll). 



We now use Lemma 4.1 again to get 

Wd2 + 1 



7(iV- 1)^ 



det 



N 

yl-Di{x) 



= det [l(di+i)(d2+i) - Tab] 
= det [l(di+i)(<i2+i) - CPi2Tij,Pi2C-i] 

iV / Ar_i 

det ld,+i 



N-l 




N-J N-l 

h a. 





N 




N 


det 


-7b 


det 






N-l 







Udi + l 

7(7V- 1)2 



det 



xl 



N 



This concludes the proof. 



(4-16) 
Q.E.D. 



Remark 4.1 Notice that this proof was based purely on an algebraic reinterpretation of the characteristic equa- 
tions for Di{x) and D_2{y) in which the same set of recursion parameters {aj(n), I3j{n)} appear. No a ssumption 
was required about any relations between these parameters, and therefore the equalities (4-7), (4-8) are really 
just identities. 
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4.2 Duality pairings 

In what follows we will derive a deeper form of duality; namely, that the linear differential equations 
satisfied by dual windows are also dual, in the sense of having the same associated spectral curves. 
This follows from the generalized Christoffel-Darboux formulae satisfied by kernels Kij,i,j = 1,2 when 
the biorthogonal polynomials and their Fourier-Laplace transforms are replaced by any solution to the 
equations of Prop. |3^ . 

Proposition 4.2 If {^^(y)}rieN and {0n(y)}neN are two arbitrary sequences of functions satisfying both 
the recursion relations under multiplication by y and the differential relations under application of dy 
(constructed as in Prop. |3.lD, then 



N . 



N 



{dy' +dy)[ * (y'), B f(y) \ = {y - y') I 3^(y'), (A)* $ (y) 



N 



N 



(4-17) 



and 



'N-l 



N 



N 



N-1 



id,' + 9.) I i ix'), A fix) \^ix~x')l Mx'), (B)* ^ ix) 



(4-18) 



Proof. Wc shall prove the equality (4-17) only; since (4-18) is proved identically. The expressions on 
either side of eq. (4-17) read (understanding the ip 's to depend on y' and the (p^s on y): 



N ~ 



idy'+dy)\^:^iy'),Mfiy)j 

(di J-1 
N-j+k 
j=l k=0 

d2 d2 

7iN-l)l^_^ E MN + l)4'N+i~liN-l) E MN-l)l^_^_ 

di J-1 d2 

+ E E E MN+miN+k)i^^_^ 

j = l fe=Oi=-l 
di j-1 d2 

- E E E PjiN+k)aiiN+k-j + l)l^^ 

j=l k=Ol=-l 



N 



(4-19) 



and 



iv' -y) i(2/'),A $ iy) 



N 



= iy' -y)\liN~l)i^ 



N 



d2 3-1 

-i-Y.Y.'^,iN+k)i 

3 = 1 k=0 



di 



N+k-j 
di 



(f'N+k 



= -,iN~l)^N-i E f^3iN+j)i^^-liN-l)l^ E P3in~^)^N-i- 



3 = -l 



3 = -l 
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d2 l — l di 

■ E E E MN+m{N+k)i^^^jN^^j . 

1=1 fc=oi=-i 



(4-20) 



The claim now is that these two expressions are the same. There are two ways of proving this. The 
first is a straightforward computation collecting all bilinear terms of the form 



Fpq{y',y) ■^±i^+pi.y')<i^N+q{y) 



(4-21) 



in the difference between (4-19) and ( |4-20 ) and proving that their coefficients vanish. The coefficient of 
Fpq with p > and q > vanish identically, as do the coefficient of Fpq with p < — 1 and q < — 1. The 
coefficient of Fpq with p > and q < —1 vanishes due to relation (2-31 ) with I — p — q — 1 and n = N + q. 
The coefficient of Fpq with p < —I and q > vanishes due to relation ( ^-32 ) with I — q — p — 1 and 
n = N + p. These cancellations are summarized in Table |^. 

The second way does not involve any computation; in fact, we can already conclude that the coefficients 
of all terms ^j) ^_^^{y')(l)N+q{y) must agree. We know that for the polynomial solutions {^„(x), (^„(?/)} 
and their Fourier-Laplace transforms {■(/' {y'), (2;')} t^i^ t'^o expressions are the same because 



N-l 



{dy. +dy){^ iy'), B $(2/) j = {dy, + Oy) (?/' ^y)^^ ' ) (^ ) 

N-l y 

= {y' - y) {dy, + dy) J2 tJy')My) = iy' - y) (my'), (a)* %\y) 



where we have used the generalized Christoffel-Darboux formulae and the identity 

[{y' -y)Ady,+dy)] = Q . 
Since the Tp^{y'ys are linearly independent and so are the (/)„(?/) 's, the functions 

Fpq{y', y) ^p_^^{y')<j)N^{y) 



(4-22) 



(4-23) 



(4-24) 



are also linearly independent. Considering eqs. ( [4-19| ) and ( |4-20 ) as linear equalities for the Fpq{y' ,y)^s, 
one concludes that the coefficient in the two equations must be equal. Q.E.D. 

Before stating the next result, we define new pairings by studying the effect of deformations on the 
kernels. By using the same Christoffel-Darboux trick one easily computes 



d 



N-l 



K 



duK 

d 



-N-l 



n=0 
Af-1 



j=l 1=1 
K J 



dv 



K 



i^n{x)±y) = J2J2 ^""(^-ot 

j=l 1 = 1 



n=0 



N- 



(4-25) 
(4-26) 



This pairing does not change if we take Fourier-Laplace transforms of 5* or (f>, so we can easily write the 
deformations of the two kernels Kn and K22- 
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Proposition 4.3 For any two sequences satisfying both tlie deformation equations and the x-recursion 
relations we have 



d 



duK 
— I 



^^-1 N . ' 

± (y'),B$(2/) 



J 3 

.j=i 1=1 



(f>N-l 



■N-l 



d N 
OUR \ N 



(^'-^)|EEt^f(^-OV^iV-/+: 

, J=l 1=1 



dv 



Q fN-l N 

N 



J j 



, j=i 1=1 



N-l 



(4-27) 
(4-28) 
(4-29) 
(4-30) 



Proof. We prove only one identity, the others being completely similar. If the two sequences consist of 
the orthogonal quasi-polynomials (and the corresponding Fourier-Laplace transforms), then the equality 
follows immediately from: 



d 



- f * (y'),B$(2/)') = -^{y' -y)y^i^iy')K{y) 
'< \N-i J duK 



n=0 
K j 



(2^'-y)|EEt^f(^"0V:^- 

, j=i 1=1 



(4-31) 



Expanding both sides by means of the recursion relations and the deformation equations in order to 
obtain linear expressions in Fp^q{y',y) in the LHS and RHS one concludes that the coefficients must be 
the same. But the same final expression relies only on the recursion relations, and hence the equality 
holds for any pair of sequences of functions ip^[y') and 4>n{y) satisfying these same recursion relations 
(by the same argument as in the proof of Prop. |4.2| ). Q.E.D. 

Theorem 4.1 If {i^^{y)}neN and {^„(y)},ieN (or {0^ (a;)}n gM and {'4'n{x)} neN) ar e arbi trary pairs of 
sequences of functions satisfying the recursion relations 1 ^^ , ( p^ ), (resp. p36|)), (|2^, th e differen- 
tial relations ( |2-25D , ( ^-21| ), (resp. ( |2-27D , ( |2-19D ) and the deformation equations ( |2-85D -( ^-88| ), then the 
bilinear expressions 



'N-l 



fNiy) 

9n{x) 



N - 

* (y),B|(y) 



N -. 



$ (a;), A^'(a;) 



AT 



(4-32) 
(4-33) 



are independent of y (resp. x) and N, and also are constant in the deformation parameters {uk,vj}. 
Proof. Using Prop. |4.2| and setting y — y' we find at once that 

|My) = o, 
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i-S- fN{y) — In does not depend on y. A similar computation shows that gN{x) is independent of x. 
Now we also compute (for, say, M < N) 



* (j/),B|(y ) 



(2/),B|>(2/ ) 



N~l 



(4-34) 



=M 



Letting y = y' we obtain /jv = /m) and similarly for ^at To prove the independence of the deformation 
parameters uj^ and vj, we use Prop. 4.3 in a similar way to the above and then again set y = y' (or 
X = x') to conclude the proof. Q.E.D. 

Theorem 4.1 allows us to conclude that we can choose fundamental systems of solutions to the pairs 
of dual differential-difference-deformation equations normalized in such a way that the pairing gives the 
identity matrix. 

Corollary 4.1 There exist two pairs of sequences of fundamental matrix solutions to the difference- 

N-l JV-l 

differential-deformation equations (|j)-(Q, ( |3-29D -(^-32D (*, * ), ($, * ) such that 



N-l N 

* ,A*) 

N 



N 

N-l N 

N 



N 



(4-35) 



We conclude with the following theorem. 
Theorem 4.2 The differential-deformation systems 



N 



dy^ = - D2{y) * = - Ua' d^, * = Vj 



N 



N 



N N N N 

N-l N-l N N-l N-lN,-i' N-l W-l W,* 

9y * = * D,{y), du^ * - * U^, 9,, * = - * V^, 



for K = 1 . . .di + \] J = l...d2 + 1, are put in duality by the matrix 

N N N N 

MD2{y)=D2{y)M , 

N NN.<S> N,<S> N 

B = B VK{y) - HAy) B 

N W.* N JVAf,* 

dv.,B,= Y,;(y)B-B Vj(y) . 



N 



(4-36) 
(4-37) 

(4-38) 
(4-39) 
(4-40) 



N N ^ 

In particular, since the matrices D^iv) and D2{y) are conjugate to each other, their spectral curves are 
the same. Similarly, the differential-deformation systems 



N N.-ii TV.* 

N N N N ' N N 

N-l N-l N N-l N-l N,<I> N-l N-l N,<S> 

9. * - * D,ix), du^ * = - * U^, d^, * = * V,, , 

N 

K — 1 . . .di + l\ J=l...fi2 + 1, are put in duality by the matrix A, 

N N N N 

ADi{x) ^Dj^{x)A , 

N NN,-^! N,iS> N 

d^jA = A Vj{x) - Yj{x)A , 
duK^^ Vj^{x)A-A IJKix) , 



(4-41) 
(4-42) 

(4-43) 
(4-44) 
(4-45) 
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N N 

and hence the spectral curves of D.i{^) and Di{x) are also the same. 



Proof. The three relations follow easily by taking a funda ment al system of solutions for the two com- 
patible differential-deformation systems and using Theorem 4.1. Q.E.D. 



This theorem together with Prop. proves that the four spectral curves 

Prop. 4.1 



det 



det 



N 

2/1-^2(2;) 



= 



I Thm. y 

N 

yl~D2{x) 



' rop. [Til 



det 



det 



xl 



N 

Di{y) 



I Thm. 4.1 



= 



(4-46) 



N 



all coincide. 



4.3 Concluding remarks 

In this work, the main results concern the compatibility of the difference- differential-deformation systems 



arising from the "folding" procedure (Proposition 3.1), and the resulting spectral duality Theorems 4.1 



and 4.2. The constancy of the bilinear pairings between solutions given by Corollary 4.1 may be viewed 



as a form of the bilinear relations for Baker functions which imply the Hirota bilinear equations for the 



associated tau function of eq. (1-2) 



Another consequence of this compatibility is the fact that the (generalized) monodromy of the covari- 

N N 

ant derivative operators ^ — Di [x) , and ^ — D2 (y) is independent of both the continuous deformation 
parameters {m_r-, vj} and the integer TV; i.e., we have a dual pair of differential operators families whose co- 
efficients satisfy differential equations in the parameters {uk, vj} and difference equations in the discrete 
parameters N that generate isomonodromic deformations. Associated to such isomonodromic deforma- 
tion equations, there is a sequence of isomonodromic tau functions in the sense of refs ||3l] , |3^ . However, 

N N 

since the highest terms of the polynomial matrices Di{x) and D2{y) have a very degenerate spectrum 
(in fact, they have rank 1), the standard definition of the isomonodromic tau function does not apply. 
To introduce a suitable definition for this situation, an analysis of the formal asymptotics at x = 00 (or 



y — 00) is required. Also, the systems of Proposition 3.1 represent in a sense, the "vacuum" isomon- 
odromic deformation systems associated with the Fredholm kernels appearing in Proposition 2.1. When 
the corresponding integral operator is supported on a union of intervals, the computation of its resol- 
vent is equivalent to a Riemann-Hilbert problem with discontinuities given across these cuts [^9|. The 
resulting "dressed" Baker functions determine isomonodromic families of covariant derivative operators 
having, in addition to polynomial parts, poles at the endpoints of the intervals, which may be viewed as 
new deformation parameters. The associated isomonodromic tau functions are given by the Fredholm 
determinants of the integral operator supported on the union of intervals. p9[ . 

The study of the formal asymptotics associated to the vacuum systems, the corresponding isomon- 
odromic tau functions and the relation of these to the spectral invariants will be developed in a subsequent 
work (1^), as will the study of the N ^ 00 asymptotics of the biorthogonal polynomials and associated 
Fredholm kernels. 
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A Appendix: Multi-matrix model 



The "multi-matrix-modcl" is a generalization of the 2-niatrix-modcl, which was introduced in the context 
of string theory and conformal field theor y p^ , , and has been extensively studied ||3^, |l^ . Our 
notations in the following mainly follow pO|. Calculations of spectral statistics in this model again 
involves biorthogonal polynomials which obey linear differential systems of finite rank. It will be shown 
in this appendix that these also satisfy an extended form of the spectral duality relations derived for the 
2-matrix case. The results will be summarized, but only a brief sketch of the proofs will be indicated. 

Consider m > 2 random hermitian N x N matrices Mi, AI2, . . . , Mm, with the measure 

?n rn—1 M 

k=l k=l k=l 

where dM^ is the standard Lebesgue measure for hermitian matrices, and the potentials Vk, k = 1 ... to, 
are polynomials of degrees d/c + 1, with coefficients 

1=1 

As in the 2-matrix case, all the correlation functions and statistical properties of the eigenvalues of 
the TO matrices can be expressed in terms of determinants involving to^ Fredholm integral kernels, which 
are constructed from an infinite sequence of biorthogonal polynomials and their integral transforms . 
In this case, the biorthogonal polynomials 

TT,,{x) ^ x'' + . . . , a„(2/) = 2;" + ... (A-3) 

are defined to be orthogonal in the fohowing sense: 

/m m — 1 

dxi dx2 ■■■ dxm-i dx,n [| e-^'=(^'=) J| e^'=^'=+i n„{xi) ai{x„,) = hj„i , (A-4) 
fc=i fc=i 

where the integral is convergent on the real axis if all the degrees + 1 are even, and the leading 
coefficients are positive. Otherwise we need to consider other integration paths in the complex plane, 
without boundaries, so that integration by parts may be done. This uniquely determines the polynomials 
TTn and an for all n. 

From 7r„, we define the following to sequences of functions {V'i.n}„^Q 00 ' ■ ■ ■ {V'ni,n}„^o 00 
i'lA^) -^7r„(a;)e~^^(") , 



■02,n(a:) := J dy V'i,n(y) e^^^ , 
^k+i,n{^) J dy ijkAv) e"" e-^''^^) for , to - 1 > > 2 , (A-5) 

and from (T„, the following to sequences of functions {'/'i,n}„^o...oo ' ■ ■ ■ {'^™."}n=o...oo 
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<jyk-iAx) y dzj 0fc,„(y) e^^ e-^-'-i^^) for m>k>2, 

<Pi.n{x) := Jdy^2Ax)e^y, (A-6) 
which are dual bases for the respective spaces they span: 

da; i^k,n{x) 4>k,i{x) = Sni ■ (A-7) 



A.l Recursion relations 

We define the semi-infinite matrices Pk and Qk for each fc = 1, . . . , m, such that ^ 

d 

Xtpk,n{^) = ^Qk(n.l)i^k,l{x) , —1pk^n{x) = ^ Pk{n,l)'^Pk.l{x) , (A- 

l I 



where it wih be shown befow that these only involve finite sums. From the pairing (A-7), and integration 
by parts, we have 

Xk(t)kAx) ='^Qk(n,l)(t'kAx) . ^0fc,n(a;) = - ^Pl(n,l)'l'kAx) ■ (^"9) 

I I 

Note that these matrices all satisfy Heisenberg relations |^ 

[Pfc,Qfc]-l. (A-10) 

Using the definitions of ^k,n and (f>k.m we find the following relationships between them. 

Pk+i ^Qk ke [1, m - 1] , -Pk = Qk+1 - V^iQk) fc e [2, m - 1] , -Pi = Q2 • (A-11) 

In particular, this implies that 

Qk-i + Qk+i - VliQk) for k e [2, m - 1] . (A-12) 

These relations are enough to ensure that all the matrices Qk and Pk are of finite band type. 

Proposition A.l The matrix Qk has bands above the principal diagonal and Sk bands below the 
principal diagonal, with 

k-l 

ri = 1 and = forfce[2,TO], (A-13) 

1=1 

s„i = 1 and Sk = '^i for fc S 1] . (A-14) 

l=k+l 

Proof: Qi multiplies the vector of polynomials {7r„(a;)}„=o...oo by x, and can therefore raise the degree at 
most by 1; i.e. Qi has at most one line above diagonal. From the same argument, since the multiplication 
by Pi + V{{Qi) takes the derivative of the vector polynomial [7r„(a;)]„=o...cx; with respect to x, it must 
lower the degree, therefore Pi has at most di = deg V( lines above diagonal and so does Q2 = —Pi- Using 



A-12 recursively, it follows that Qk has at most lines above diagonal. Repeating the argument with 



the polynomials {o-„}, we see that Q/. has at most Sk lines above the diagonal. Q.E.D. 
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Denoting 

ak,i{n) := Qk(n.n+i) 
the recursion relations may be written componentwise as 



(A-15) 



vipk,n{x) ^ ^ ak,iin)'ipk,n+iix) , x4>k,n(x) ^ ^ ak,i{n - l)4>k,n-l{x) , (A-16) 



l = — Sk 



which impUes 
d 



+r2 



dx 

and 
d 



d 

^ a2,ii^)'^i,n+i{x) , -^tl^k,n{x) ^ ^ ak-i,i{n)ij)k,n+i{x) , k ^ [2,m] , (A-17) 

l = — S2 



+r2 



d 



+rk- 



— (/)i^„(a:) = ^ a2,i{n~l)(l)i,n^i(x) , j^(j)k,n{x) = - ^ afc_i,i(n - Z)(/)fc_„_;(x) , /c € [2, m] . 



Note that 



l = ~S2 



dx 



ai,i(n) = a,n-i{n + 1) = 7('^) = 



hn+l 



(A-18) 
(A-19) 



A. 2 Folding 

Again, it is possible to "fold" these recursion relations to form finite rank linear differential systems with 
polynomial coefficients. For each fc, define the following windows of size + Sk- 



N 



k,N-Sk 



^k = ( 4'k,N-rk 
N 



'^k,N+Sk 



.-l) 



\1pk,N+rk-l , 

Shift operators: For each k, define "ladder" matrices of size rk + Sk, 



ak{x) 

N 



ak{x) 

N 







V 











Ck.-sJN) 



x-ak,aiN) 
ak,r,.{N-rk) 



ak.-sJN+Sk) 
ak,r^{N-rk) 



■•■ 

1 




which implement the shiftsR in N. 



N N N+1 



^kClkix) = $fc 

AT N N-1 



^In the m = 2 case, we had a = ai, a = iii, b = (aj) b = (aj) ^ ■ 



(A-20) 



(A-21) 



(A-22) 
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Differential systems: The differential systems^ satisfied by the vectors 'i'kNi^) ^^'^ ^kwi^) ^'I'e: 
g N N N r) !y 

— ^kix)=Dk{x)^k{x), j-<^k{x)^~^k (x) Dk (x) (A-23) 

ox ox N N 



N 



N 

where Dk{x) and Dk{x) are matrices of size + Sfc, and with polynomial coefficients of degree at most 
dk in X. 

N 

The matrices Dk{x) are given by: 

jY r2 I S2 l-l / \ -1 

£'1(0;) = - a 2,0 - a 2,; TT ai (a;) - a 2,-i TT I oi (a;) 

^ fcAT >. - kN T-i- X - kN Tf I 

Dk{x) ^ afc-i.o+2^ "fe-i.iJ_J_ flfc (2;) + 2^ "fc 

1=1 0=1 /=i j=o V^-'+J' 

for m > fc > 2 , (A-24) 

where 

fcAT 

a J- i := diag (a^- i(iV - Sfc), . . . , a^- i(iV + rfc - 1)) . (A-25) 



Similarly: 



N ru-i I l-l \ Sk-i / I 

Di{x) = - a 2,0 - 51 n ^1 " - 51 n ( "1 .(^) 



a 

2,-; 



k-l,-l 



Dk{x) = a k-ifl + TT (^) " fc-i.' + TT ( (2^) 

for m>k>2 (A-26) 



where 



a J- i := diag (a^- " I " '^fc), • • • , aj,;(iV ~ I + Sk - 1)) . (A-27) 



Christoffel Darboux matrices: 

Consider the kernel 

Af-l 

K ix,y) = V 4'kAy)^kAx) ■ (A-28) 

Nk,k 

The generalization of the Christoffel-Darboux theorem for this kernel reads 

rk I Sfc I 

{x~y)K {x,y) = ^^ak,i{N ~ j)^k,N-j+i<l)k,N~j -^^ak-i{N - j + l)i}k,N-j<l>k,N-j+i 
^ 1=1 j=i 1=1 j=i 

= $fc(x)Afc^'fc(y) . (A-29) 

N N N 



® Notice that the notations are changed for m = 2. D2 is now what we called —Dj^ and D2 is what we previously called 
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where is the (r^ + Sk) x (rfe + Sfe) matrix: 

N 

/ 








V 






-ak,-i{N) 



ak,rk{N-rk) 



aM(^-l) 






\ 



afe,r,(A^-l) 



/ 



-ak^+s,-i{N + Sk ~ 1) 

There is also a "differential Christoffel-Darboux theorem", for k > 1: 

Tk-l I 

{dx + dy)K {x, y)= ak-i,i{N - j)ijjk,N-j+i(t>k,N-j 

Sk-l I 

- X] 5Z - 3 + l)ll)k,N-i(l>k,N-j+l 

1=1 j=l 

= $fc(x)Afe_i#fe(y) , 

N N N 



where 



^k = {'4'k,N-Sk-i ■■■ V'fe,iV+rfe_i-l ) , ^k = {4'k,N-rk-i ■■■ 0fe,iV+Sfc_i-l ) 

AT JV 



(A-30) 



(A-31) 



(A-32) 



and for A; = 1, 



where 



r-2 I 

{dx + dy)K {x,y) = -^^a2,i{N - j)tpi^N-j+i(pi,N-j 

S2 I 

,N-j4>l,N-j+l 

1=1 j=l 

= -$i(x)A2*i(y) , 

AT JV AT 



*1 = ('^i'l,JV-S2 ••• ■^i'l.JV+rs-l )* ) ^1 = ( <?^l,JV-r2 ••• <l>l,N+S2-l) ■ 



JV 



JV 



(A-33) 
(A-34) 



A.3 Duality 

JV JV 

All the systems Dk{x) and D^{x) have the same spectral curve. This result follows again in two steps. 
• For each k, we have: 



JV 



N 



Dk{x) Ak = Afe Dk{x) 

JV JV 



(A-35) 
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which imphes that 



N 



N 



det yl - Dk{x) = det yl - Dk{x) 



The relationship between spectral curves for different k is: 

det (x\ - Dl+i{y) \ cx det Ul - V;,{x)l + Dk{x] 



and 



N \ f N 

det ( xl - D2{y) oc det yl + Di{x] 



for fc > 1 



(A-36) 

(A-37) 
(A-38) 



Proof: We prove A-35| using the same method as for m = 2. Here is a sketch of the proof for fc > 1. 
Let ^kN{x) and ^kNiu) be any solutions of the differential systems 



— *fe(x) ^Dkix)'^kix) , 7r*^-(y) = -*fc(2/)^fc(y) 

ox N N Oy N N 



(A-39) 



We construct the functions ipk,n{x) with N — Sk — Sfe-i < n < N + rk + r^-i — 1 and 4>k,n{y) with 
N — Tk ~ Tk-i < n < N + Sk + Sfc-i — 1 by recursively applying the shift operators (which are again 
compatible with the differential systems). This gives 



{d^ + dy) ^k{y) Afc ^fc(x) ^{x-y) ^k{y) Afc_i #fc(a;) 



N N N 



N 



N N 



(A-40) 



This equality holds term by term, since the coefficients for each monomial of type ipk,nix)(j)k,i{y) are the 
same as when ipk.nix) = '4'k,nix) and 4'k,i{y) = 4'k,i{y)i which are linearly independent functions of x and 
y. By taking x = y one has, for any ^kNiu) and ^'^jvl^;), 



N 



N 



$fc(a;) Dk{x) Ak - Ak Dk{x) = . 

N \ N N IN 



(A-41) 



Since this holds for any basis of solutions, the factor in brackets (. . .) must vanish. 



Q.E.D. 



Here is a sketch of the proof of ( A-37 ); it is very similar to the proof of Prop. iA. First, we show 
how to prove that 

det (yl + Di{x)) oc det {xl - D2{y)) , (A-42) 

the other cases being similar. 
First, notice that: 



det {yl + Di{x)) 



X det I 1 — fli 

JV-S2 



det (-'0^2,-82) 



det I ai . . . fli 



IN 

a 2,-s 



fli . . . Oi 

N-1 N-S2 + I 



Ea-i - 
AT \. 1 



l = -S2 + l 



N-S2 



IN 

a 2,1 

N 

a 2,-s; 



fli ... Oi 

N+l-1 N-S2 + 1 



(A-43) 
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Using lemma 4.1, the last determinant can be written as the determinant of a block matrix Ti of size 
in + si) X (r2 +S2). 

det (yl + Di{x)) = a det (1 - Ti) , a = const. , (A-44) 

where 


































ai 

N-S2 + 1 


IN 

a 2,r2 


IN 

2,r2-l 


IN 

y- "2,0 


IN 

a 2,-s2+i 


" 2, — S2 




" 2,-S2 





(A-45) 



On the other hand, by the same argument, we have 

det (xl — D2{y)) = C2 det (1 — T2) , C2 = const. 

where 



(A-46) 



To := 






02 

TV+r-i-l 




























02 

TV-si + l 


2N 


27V 

a i.n-i 


27V 

X ~ a I Q 


2Af 


°2 2W 




■■■a2N-s^ 2Ar 


°2 2Ar 



(A-47) 



It is easy to see that Ti and T2 are equal up to permutations of rows and of columns, and therefore they 
have the same determinant. 



The other equalities with fc > 1, 



N \ / N 

det xl - Dk+iiy) cx det Ul - V^ix)Id + Dkix] 



for fc > 1 , 



(A-48) 



are proved by the same method and by induction on fc. We define the sequence of functions Xj{x,y), 
1 < j < + 1; such that 



Xj-i = V-{xj) - Xj+i 2 < j <k . 



(A-49) 



Xk = X , Xk+i = y 

We then prove by induction on j that 

det {xj^i — Dj{xj)) cx det (xj — _Dj+i(a;j+i)) 2 < j < k . (A-50) 

Each step of the induction is similar to the method described above for fc = 1. This completes the proof 

Q.E.D. 



of A-37 



It can also be proven that all these systems are compatible with the shifts and deformations. It follows 



that if $fc(x) and *fc(a;) denote fundamental solution matrices for the systems (A-23), it is possible to 

N N 

choose their normalizations such that 



$fe(a;)Afe^'fe(.T) = 1 . 

N N N 



(A-51) 
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This may again be viewed as a form of the biUnear identities that aUow us to deduce bihnear equations 
for r-functions. 
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CoefT. of 


Eq. (4-19) 


Eq. (4-20) 




7"(iV-l) + 
1 a,fiV + 7 - l)0i(N + j - 1) 


7^(iV- l)-f 
y'^ , ai(N + j - DBAN + 7 - 1) 




-7'(iv-i)-E-iia.(mW 


-i\N-i)-j:%aANmN) 


p>0,q>0 
p + q>l 


- ^ a,(Af + g)/3p-,+i(A^ + p) 


- aiiN + q)Pp^^+i{N + p) 
l=g+i 


P + 9 > 3 


^Qfe+p(7V + fc)/3fc+,(iV + fc) 

fc=0 


J2ak+piN + k)Pk+g{N + k) 

k = 




7(7V - l)ao(iV) - 7(iV - l)ao(iV - 1) 


-^a,(7V)/?,_i(iV-l) 




^a,_i(7V + j-l)/3,(iV+j-l) 
j=i 

~^a,_i(7V-l)/3,(iV) 

3 = 1 


7(iV - l)f3o{N) - 7(7V - 1)Pq{N - 1) 


q>2 


^a,(iV + j)/3.+,(iV + j) 
j=o 

- ^ /3,+,(iV)a,(iV-g) 
i=-i 


-7(iV-l)/3,_i(iV-l) 


p > 1 


^ a,(iV+p + j)ft+P+i(^' + P + j) 
j=-i 

-^a,(iV-l)/3,+p+i(iV + p) 

J=0 


7(iV-l)/3p(iV + p) 


9> 1 


7(iV- l)a,(A/- + q) 


^a,(iV + j-l)/3,„,_i(iV + j-l) 
-^a,+,+ i(7V + Q)/3,(7V-l) 

j=0 


In-p^^ 
p>2 


-7(iV-l)ap_i(iV-l) 


^a,+p(iV + i)/3,(iV + i) 

- ^ a,.+,(iV)/3,(7V-p) 
i=-i 


±N+p^^-^ 

q>'2, p>o 


^ (iV + J + p)/3,-+p+, (iV + j + p) 

- Pj+P+<,{N + p)aj{N-q) 
j=-i 





In-p^^+'' 

g > 0, p>2 





^ a,+,(iV + j)/3,_,(iV + j) 

- aj+p+,{N + q)l3j{N-p) 
j=-i 



We have defined aj := if j ^ [—1, . . . , c?2] and Pj = if j ^ [—1, . . . , rfi] 



Table 1: Comparison of coefficients in eq.(4-19) and eq.(4-20) 
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